% test_spectrum2D
s1 = lpcvd_2DSurface(filename);
sp1 = s1.cal_spectrum();

s2 = s1.aggregate(10);
sp2_a = s1.cal_spectrum();
sp2_b = sp1.aggregate(10);

fprintf('From definition: %f, From mode: %f, %f\n',s2.r2,sp2_a.r2,sp2_b.r2)

